Characterizing the entanglement of symmetric many-particle spin-1/2 systems 
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Analyzing the properties of entanglement in many-particle spin-1/2 systems is generally difficult 
because the system's Hilbert space grows exponentially with the number of constituent particles, N. 
Fortunately, it is still possible to investigate many-particle entanglement when the state of the system 
possesses sufficient symmetry. In this paper, we present a practical method for efficiently computing 
various bipartite entanglement measures for states in the symmetric subspace and perform these 
calculations for N ~ 10 3 . By considering all possible bipartite splits, we construct a picture of 
the multiscale entanglement in large symmetric systems. In particular, we characterize dynamically 
generated spin-squeezed states by comparing them to known reference states (e.g., GHZ and Dicke 
states) and new families of states with near-maximal bipartite entropy. We quantify the trade-off 
between the degree of entanglement and its robustness to particle loss, emphasizing that substantial 
entanglement need not be fragile. 



I. INTRODUCTION 

The structure of entanglement within multipartite 
quantum systems is a deep subject that has only be- 
gun to be explored. Because an ensemble's Hilbert space 
grows exponentially with the number of particles that 
comprise it, the number of distinct ways these particles 
can become entangled and the number of reference states 
needed to represent the various entanglement structures 
are immense While exponential scaling in complex- 
ity is the reason multipartite entanglement is so rich, it 
is also the reason the subject is so daunting. 

Nonetheless, there is motivation for characterizing en- 
tanglement in many-particle systems such as atomic spin 
ensembles because of recent experimental progress in cre- 
ating and manipulating macroscopic quantum states. In 
particular, highly correlated atomic ensembles, such as 
spin-squeezed states Q, have been demonstrated SEJQ 
and promise advances in atomic interferometry |fj and 
quantum communication 0. They also provide experi- 
mentally accessible systems for studying quantum mea- 
surement, feedback, and control @. 

Spin-squeezing is intimately linked to the structure of 
the entan glem ent between individual members of the en- 
semble |9lll0f. However, without a complete microscopic 
picture of this entanglement, only limited claims about 
the structure of these correlated states can be made. In 
certain cases, an iV-spin system can be characterized as 
either entangled or separable by measuring (computing) 
expectation values of total ensemble operators jllL Il2| . 
For example, if the spin-squeezing parameter for a N- 
spin state (with polarization along z and minimal vari- 
ance along x) is less than unity 
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then the state is guaranteed to be inseparable. However, 
at this level, limited detailed information about internal 
entanglement and its robustness to particle loss 0, 0] , 
or other types of decoherence 0, 0, is available. 
In other words, entanglement tests using total ensem- 
ble operators cannot completely characterize the trade- 
off between the available entanglement resources and the 
state's fragility. 

Unlike several multipartite techniques that have been 
introduced (e.g., the N -tangle [jjj ). we approach the 
problem of analyzing iV-particle entanglement using only 
bipartite measures. Although a single bipartite split of 
a large system is rarely sufficient to characterize multi- 
particle entanglement, combining the results from many 
different splits of the system paints a reconstructed pic- 
ture of the many-particle entanglement. Furthermore, 
by repeating the analysis after removing particles from 
the system, it possible to systematically characterize the 
entanglement across all size scales and its robustness to 
particle loss. Our approach has the advantage that it 
relies upon well-defined entanglement measures that are 
both computable and physically motivated. 

Since substantial insight, and often a good starting 
point for more rigorous analysis, can be gained from nu- 
merical simulations, an efficient means of calculating en- 
tanglement measures is desirable. Section IIIII develops 
the necessary machinery for calculating these measures 
in the symmetric subspace — the set of those TV-particle 
pure states that remain unchanged by permutations of 
individual particles 0, 0, 0] . The main result of this 
section is that it is possible to perform partial transposes, 
partial traces, and Schmidt decompositions of symmetric 
states without resorting to an exponentially large repre- 
sentation of the system. 

In Section ITVl we characterize microscopic entangle- 
ment and its robustness to particle loss for several repre- 
sentative symmetric states, including the GHZ and Dicke 
(e.g., W) states. Here, the advantage of exploiting sym- 
metry is clear; we perform entanglement calculations for 
systems with N ~ 10 3 . These numerical results allow 
us to speculate on the large N asymptotic scaling of the 
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above entanglement measures. In some cases, particu- 
larly for the entanglement of formation and the reduced 
entropy, we analytically verify the observed scaling. We 
also introduce a family of states that provides insight into 
the scaling of bipartite entanglement in symmetric states 
for large N. 

With the context provided by the reference states and 
the boundaries of allowed entanglement structures, we 
can better understand the entanglement generation abil- 
ities of certain dynamical processes. SectionlVlfocuses on 
the entanglement produced by spin squeezing Hamiltoni- 
ans. We illustrate the intuitive and generic effect that 
small scale correlations peak before (and transform into) 
larger scale correlations. Again, the ability to simulate 
systems with N 3> 1 permits us to determine asymptotic 
behavior, both for large numbers of particles and for long 
times. 

A point we stress is that significantly entangled states 
need not be fragile. Robustness is critically important 
in experiments where the system constantly exchanges 
atoms with the surrounding environment. Moreover, we 
show that spin squeezed states provide a reasonable com- 
promise in this trade-off; they are highly entangled, yet 
particularly robust. 

II. ENTANGLEMENT MEASURES 

In this section, we review several common entangle- 
ment measures as motivation for the symmetric state 
techniques that are developed in Section IlIII In addition 
to recognizing the specific operations necessary for com- 
puting these entanglement quantities, we also describe 
their strengths, weaknesses, and, where possible, physi- 
cal motivation. 

We begin by reviewing the commonly accepted set 
of properties that all measures of entanglement should 
share. For a general density matrix, p, which can be di- 
vided into two or more sub-systems, the quantity Ex (p) 
(the label X is used to denote a generic measure) qualifies 
as an entanglement monotone if it satisfies the conditions 

(CI) E x (p) > 0; E x {p) = if p is separable; 
£x(Bell State) = 1. 

(C2) Local operations, classical communication and 
post-selection (LOCC) do not increase Ex(p) on 
average. For example, with any state, p, and 
partition, {A, B}, local unitary transformations, 
U = Ua <8> Ub , do not affect E x {p) ■ 

(C3) Entanglement is convex under discarding informa- 
tion T,iPi E x(pi) > Ex(E,iPiPi)- 

We define the generalized Bell States as, 

I**) = (|1a0u>±|0a1b»/V2 (2) 
= (|U1b) ± |(U0b))/V2 (3) 



for a partition {A, B} |^. If the subsystem A has more 
than one spin, 1a is interpreted as li • • • 1jv a an d simi- 
larly for 1b, 0a, and b ■ 

A. Entropy of Entanglement 

Given a pure state, and a partition for the system, 
{A, B}, the entropy of entanglement is defined as, 

E(\*),{A,B}) = S(pa)=S( P b) (4) 

where the von Neumann entropy is S(p) = — Tr(plog 2 p) 
and pa = TrB(\^)(^\)- Any entropy that results from 
performing a partial trace on the system must be a con- 
sequence of initial entanglement provided that the initial 
state is pure. For product states, |^) = \^}a <8> \^)b, 
the entropy is zero since the single eigenvalue for each 
of the pure states pa and ps is one. The maximum 
entropy of entanglement given a partition with dimen- 
sions, dim(yl) = dA and dim(_B) — ds, with dA < ds, is 
log 2 (rfy!t). A state that achieves this maximum is, 

I*) = |0)a®|0)b + |1)a®|1)b + "- 

+\ d A-l)A®\dA-l)B (5) 

The entropy of entanglement has the attractive feature 
that it is straightforward to compute; it requires only 
performing a partial trace, pa = Tr b(p), then computing 
eigenvalues of the result. The drawback of the entropy 
is that it only qualifies as an entanglement monotone for 
initially pure states. 

B. Entanglement of Formation 

The entanglement of formation [23] is defined as 

E F (p,{A,B})= min J2 Pi E(\iPi), {A, B}) (6) 

{pi,if>i} 

where the {Pi,ipi\ satisfy the condition that p = 
TliiPi\ 1 Pi){' t Pi\- This quantity is difficult to compute for 
mixed states but reduces to the entropy of entanglement 
for pure states. 

In the special case of a mixed state of two spin- 1/2 par- 
ticles, the entanglement of formation can be computed 
from the two-particle concurrence, C(p) [23112^. There- 
fore, it is generally possible to compute the entanglement 
of formation between two spins removed from an 

TV-spin state \^f) . The entanglement of formation for such 
a reduced system is a strong measure of the robustness 
of that state's entanglement to particle loss. Explicitly, 
for the two particle state p — Tr^ij |^')(^ r | 



E F (p, = h(-[l + Vl-CW 2 ]) (7) 

where h(x) — — a;log 2 (a;) — (1 — x) log 2 (l — x) and 

C(p) = max(0, VM- y/% - y/% - y/xl) (8) 
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in which Ai, A4 are the eigenvalues of p{o~ y ®o~ y )p* (o~ y ® 
a y ) in decreasing order and <r y is a Pauli spin matrix. 

C. Distillable Entanglement and Negativity 

Given a mixed state, p, and a partition, {A,B}, the 
entanglement of distillation is defined as, 

TCI 

E D (p,{A,B})= lim - (9) 

n — >oo n 

where m is the number of Bell states that can be distilled 
from n copies of p via an optimal purification protocol 
with LOCC l2r| . For simplicity, we consider only 
the symmetric Bell state |$ + ) of Eq. as the output 
of the distillation process throughout this paper. This 
state is also known as an EPR pair, a GHZ state, or 
an TV-particle cat state. The distillable entanglement is 
effectively a conversion efficiency; however, since the pu- 
rification protocol allows auxiliary separable states to be 
introduced into the original system, it is possible, on av- 
erage, to extract more than one EPR pair from an ini- 
tially entangled state. The distillable entanglement for 
an EPR pair is one by definition. 

The advantages of the distillable entanglement are that 
it is a monotone for mixed initial states and that it quan- 
tifies entanglement as a practical resource. In this sense, 
the distillable entanglement has a direct physical inter- 
pretation. Unfortunately, it is extremely difficult to com- 
pute unless the initial state is pure, in which case it re- 
duces to the entropy of entanglement. The entanglement 
of formation is an upper bound on the distillable entan- 
glement (i.e., one cannot extract more EPR pairs than 
the number used to form the state). 

There exists another entanglement monotone, the log- 
arithmic negativity, which, like the entanglement of for- 
mation, provides an upper bound on the distillable en- 
tanglement but is also computable for mixed states |27| - 
The logarithmic negativity is defined as 

E x {p, {A, B}) = log 2 (2AA(p, {A, B}) + 1) (10) 

where Af(p, {A, B}) is the negativity of the state, p. The 
negativity is defined as the absolute sum of the negative 
eigenvalues of the partial transpose with respect to A, 
P Ta . So 

i 

where Ai are all of the eigenvalues. 

The logarithmic negativity can be directly computed 
from the partial transpose. However, both the logarith- 
mic negativity and the distillable entanglement are zero 
for those entangled states with positive partial transposes 
(PPT). PPT enta ngle d states and perhaps some other en- 
tangled states [23, H3 have zero distillable entanglement 
|30|. These states are known as bound entangled states. 



As with all monotones, the negativity may also dis- 
agree with other monotones, like the entanglement of 
formation, on which state of two is more entangled |22j . 
This ordering problem is a caveat which qualifies many 
statements about entanglement, and is a reflection of the 
fact that any given entanglement measure refers only to 
its own limited physical context. 

D. Schmidt Decomposition 

For a given partition, {A, B}, of the full ensemble's 
Hilbert space, it is possible to decompose the state as 

i*) = EE c ^>^ ( 12 ) 

iGA j£B 

where the kets, {|i)A)|j)s}i provide complete bases for 
A and B, respectively. For separable pure states, the 
matrix, c, which is not necessarily square, is rank one, 
R(c) = 1. States where R(c) > 1 are entangled because 
they cannot be expressed as a single tensor product. 

Generally, the Schmidt basis is taken to be diagonal 
in A. It can be found from the matrix elements, Cy , by 
performing a singular value decomposition of c, 

c = UAV 1 (13) 

where A is diagonal and the rows of U provide the 
Schmidt basis [32|. There are r = R(c) nonzero elements, 
Ai, . . . , A r , along the diagonal of A. 

Several bipartite entanglement monotones can be de- 
fined as functions of the Schmidt coefficients [3^, l34| . 
however we present this formalism only because the 
Schmidt decomposition provides an efficient procedure 
for computing the entropy of entanglement. Starting 
with a pure state, the reduced entropy for the partition, 
{A, B}, is given by, 

r 

£(|*),{A,5}) = -]TA2log 2 (A2) (w) 
i=i 

where the Ai are the singular values from Eq. I|13|) . 

III. SYMMETRIC STATES 

The previous section provided motivation for comput- 
ing partial traces, partial transposes, and Schmidt de- 
compositions. However, for arbitrary TV-particle spin-1/2 
ensembles, these operations are exponentially difficult to 
compute because a general state of the ensemble resides 
in the space <Cf N and the dimensions of the density ma- 
trix scale as 2 N x 2 N . Computational investigation of 
arbitrary ensemble entanglement is therefore impractical 
for all but the smallest values of N. 

Fortunately, a large number of experimentally relevant 
states possess symmetry under particle exchange and this 
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property allows us to significantly reduce the computa- 
tional complexity. A large class of iV-particle states are 
invariant to symmetry transformations of the permuta- 
tion group, 



Uijp N Il 



Pn, 



V IL 



(15) 



where the Ily are operators that exchange particles i and 
j within the ensemble. This is the most general class of 
states that are exchange invariant; however, it is also 
possible to further restrict the space of accessible states 
to those that are symmetric with respect to single-sided 
permutations, 



IlyPjv = Pn, V Hi 



(16) 



of the individual spins. This symmetry further constrains 
the diagonal terms of the density matrix. For the exam- 
ple of a two spin system, single-sided symmetry requires 
(01|p|01) = (10|p|01), while the more general double- 
sided symmetry does not. 

The states, \m,N), that respect this single-sided per- 
mutation symmetry compose the symmetric subspace, 
§at. The ket, \m,N), is defined as the unnormalized N- 
particle symmetric state with m excitations (spins up), 

\m,N) =^Pi(|l l5 l 2 ,...,l m ,0 m +i ) ...,0 JV » (17) 



where {Pi} is the set of all 



N 



distinct permutations 



of the spins. Although each \m, N) is an element of (D| 
the permutation symmetry enables it to be expressed 
as an element, |to), of a space, Sjv, that scales linearly, 
rather than exponentially, with the number of particles. 
In short, all states in §n can be represented in (Djv+i- 

The symmetric subspace therefore provides a conve- 
nient, albeit idealized, computationally accessible class 
of spin states relevant to many experimental situations 
(such as spin squeezing) . Completely symmetric systems 
are experimentally interesting largely because it is often 
easier to non-selectively address an entire ensemble of 
particles rather than individually address each member. 
Of course, there are still technical challenges in preserv- 
ing perfect symmetry among the particles in an ensem- 
ble, such as maintaining the uniformity of magnetic and 
optical fields. Still, for a system of many particles, sym- 
metrically manipulating the ensemble generally requires 
fewer resources than addressing individual members. 

It is therefore attractive to consider computing various 
measures of entanglement and simulating the system's 
dynamics using symmetric states. However, analyzing 
entanglement requires at least the operations of partial 
traces and partial transposes. In order for these opera- 
tions to be practical for large N, it is essential to compute 
them in an efficient manner, i.e., without having to work 
with representations of states in the full space, <Cf N . 

In this section we derive relationships that allow us 
to work with arbitrary bipartite splits of the symmetric 



subspace. The ability to express a symmetric state in 
terms of tensor products of smaller symmetric states is 
a critical prerequisite for efficiently computing bipartite 
entanglement measures. In Section IIII Al we derive the 
necessary expressions for expressing symmetric states in 
reduced dimensional bases. These results lead to the op- 
erations of partial traces, partial transposes and Schmidt 
decompositions on symmetric states. In all of these cases, 
it is possible to manipulate symmetric states with at 
worst polynomial scaling of the required computational 
resources. 



Symmetric Change of Basis and Decomposition 
Operators 



When working with the symmetric subspace, it is nec- 
essary to convert between the large, (Df N , and small, 
Cjv-fi, basis representations of the state. In order to 
provide a systematic means for changing bases, it is con- 



venient to define a symmetry operator, S_ 



N 



€ 



»N 



(Djv+i, whose action on the density operator in the 2 N 
dimensional basis, 



Pn — SnPnS n 



(18) 



projects the state into Sjv expressed in an (N +1) di- 
mensional basis. We have adopted the notation that jojv 
is the symmetric density matrix represented in (Djv+i- 

Sn is an (N + 1) x 2^ dimensional matrix that can be 
expressed as, 



JY 



Sn = ^ c N,m\m)(m,N\ 



m=0 



where the coefficients are given by, 



C 



N.m 



m\(N - to)! 



(19) 



(20) 



and CN,m\m, N) is the normalized version of Eq. 117fl . 
The state \rh) is physically the same as the 2 N dimen- 
sional state \m,N) (both have to spins up), except that 
|m) is normalized and expressed in the (N + 1) dimen- 
sional basis, 



(m\n) = 5, 
S N C N . m \m,N) 



TO). 



(21) 
(22) 



It should be noted that Sn is not a permutation oper- 
ator, but rather a projector. Therefore, it is only appro- 
priate to operate on symmetric states with Sn as, 



SnS n — Isym 

S N S N + 1 full 



(23) 
(24) 



where l sy m is the identity in the (N + 1) dimensional 
symmetric basis and lf u n is the identity in the 2 N dimen- 
sional full basis. Consequently, S\{SnPnS n Sn = Pn 
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only if pm is symmetric. Acting on a non-symmetric 
state with Sn and S N results in a loss of information, as 
the non-symmetric components of that state are lost in 
the projection onto §at. 

For the purpose of making a bipartite split, {A,B}, 
the essential property of the symmetric subspace is that 
it can be expressed as a tensor product of smaller sym- 
metric spaces. However, the tensor product of arbitrary 
symmetric states is not necessarily symmetric, 



-'.V 



C S 



N-k ' 



(25) 



where the partition {A, B} has been denoted by the num- 
ber of spins in each subsystem, {N — k, k}. §Ar_fc <g> is 
larger than Sjy. The structure of valid symmetric prod- 
ucts is given by the relation [Tflj . 



k 

|m, N) = ^2\m - p, N - k) ® \p, k) 

p=0 



(26) 



in terms of constituent symmetric states expressed in the 
large basis. 

Eqs. H25f) and l|26|l raise the point that the N particle 
symmetric space, Sjv, is smaller than the product space, 
§>N-k <8> §>fc- Therefore, the entanglement of states in Sat 
will generally be more restricted than those in the tensor 
product space. While, it is straightforward to identify the 
maximal entanglement bounds for states in, §w-fc <S> Sfc, 
the same is not true for §n- Therefore, it is convenient to 
use the product space entanglement bounds as an upper 
limit, albeit an overestimate, for the scaling of states in 

In order to exploit the tensor product structure in Eq. 
(|26|l . motivated by our desire to consider bipartite en- 
tanglement measures, it is beneficial to construct a new 
symmetry operator, TN-k,k, that maps symmetric states 
into the tensor product structure imposed by the par- 
tition, {N — k,k}. In order to be useful for computa- 
tions, both §Ar_fe and must be expressed in their re- 
spective small bases. That is, we require the mapping 
Ttf-k,k '■ Cat + i — > CAr_fc + i (g> Cfc+i. 

Constructing the operator, TN-k,k> can be accom- 
plished by decomposing Sn according to Eq. I|26[) . 



N 



Sn = ^2c N , q \q) 

q=0 



J2(q-p,N-k\(g>(p,k\ 

.p=o 



(27) 



and then operating on the expanded S N with both 5jv-fc 
and Sk, 



T N - 



N min(q.k) 



back to the large basis, extracting the tensor product 
structure, and then reducing the dimensions of the sub- 
systems down to their respective small bases. 



B. Partial Traces in the Symmetric Subspace 



In this section we derive an expression for, 

pN-k = Tr fe [pAr] 



(29) 



that avoids expressing any of the density matrices (in any 
intermediate step) in their large bases. The structure 
of the operator, T^-k,k, immediately indicates that this 
is possible since symmetric states can be expressed as 
tensor products of lower-dimensional symmetric states. 
Once the symmetric system has been partitioned, the 
partial trace is immediate. 

Although the operator, TN-k,k, can be directly applied 
to Pat, this approach condenses several intermediate steps 
that might be useful when performing calculations. In- 
stead, we first convert pn back to the large basis, 



Pn 
Pn 



S n pnSn 



N 

E 

m.n— 



(m,N\p N \n,N)\m,N)(n,N\ 

f~i-2 (~i-2 



(30) 
(31) 



and then partition the symmetric states, \m, N) and 
\n,N), using Eq. (|26[1 with k = 1. Taking the partial 
trace of the resulting expression leads to an N — 1 parti- 
cle symmetric state in the large basis, 

N 

Tri[pAr] = C N>m C Nfn {m,N\p N \n,N) 

m,n— 1 

x [\m,N- l)(n,N- 1| 

+ \m-l,N-l){n-l,N-l\] (32) 

which can be changed to the small basis using the oper- 
ators, Sjv-i and S N _ X , 

(a\p N -i\b) = C N 1 _ la C N 1 _ lb (a\p N \b)C Nia C N ^ 

+ {a^l\p N \b+l)C Nta+1 C Ntb+1 (33) 

By induction, it can be shown that the result of tracing 
k particles out of the system is, 



k 

ST/ 



(a\p N - k \b) = ^(a + j\p N \b+])Cj* °»*U C »™ 

r-T J L^N-k ,a^N-k, 

3=0 



c 



k,k 



N,q 



q=0 



p=0 



CN — k,q—pCk,p 

(28) 

to produce the necessary mapping. Here, \m)N-k G 
C-N^k+i denotes symmetric states in the subsystem, A, 
and the \n)k € C^+i are symmetric states in B. Eq. I|28|) 
has the interpretation of taking an |m) € §at, changing 



Q-p) N-k ® \p)k(q\ which resides within CAr_ fe+ i . 



b 

(34) 



C. Partial Transposes in the Symmetric Subspace 

The structure of Tw-k,k demonstrates that the partial 
transpose of symmetric states with respect to k particles, 
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fix, resides in the space §N-k ® , but not §at. There- 
fore, the partial transpose involves matrices that belong 
to Ck+i ® Civ-/s+i ; an d computing scales quadrati- 
cally in N. 

As with the partial trace, the operator T^-k,k can 
be directly employed to obtain the partial transpose; 
however, this approach hides several useful intermediate 
steps. Instead, a more explicit derivation involves trans- 
forming pm back to the big basis and employing Eq. (|26|l . 
The partial transpose, 



Pn 



N k 

E E C Nim C Nin {fh\p N \n) [ (35) 

m,n— p,q— 

| rn - p, -/V - fc)(n - g, N - k\ ® |g, fc)(p, A;|] 
can be expressed as a tensor product, 

E^v-.^r (36) 



Pn 



p,q=0 



where 



AP-.q 

A N-k 



and 



N 
m,n— 

x |m — p, TV — ft) (n — g, TV — ft| 



^'' = |g,fc)<p,A| 



(37) 



(38) 



Returning to the small basis is accomplished by evaluat- 



ing, i*^ = 5^_ fe A^_,5j,_ fc and = S k B™Sl to 



give. 



(S|i p /„ fc |6) = 



C N ,p+aC N , 



q+b 



CN~k,aCN-k,b 

Ck,lCk,d$q,c5p,d 



where, 



Pn 



= E A i 

p,q=0 



N-k 



a+p\p N \b + q) (39) 
(40) 



(41) 



shows that the dimension of p^f is, in fact, (k + 1) X (TV — 
fc + 1). 



D. Schmidt Decomposition of the Symmetric 
Subspace 

It is quite simple to perform the Schmidt decomposi- 
tion, Eq. 1|12|) . of a symmetric state in Sjy, into the space 
§v-fe <E> Sfe. The coefficients, c, in Eq. (|12|l for the states 
|m) follow directly from applying the operator T^-k.k to 
I to), resulting in the expression, 



N-k k 



Tjv_fc,fe|TO) = E E s. 



c 



N,: 



i=0 j=0 



K)iV-fe ® |j)fe 



(42) 



For the states, |m), the Schmidt matrix, c,is sparse and 
the singular value decomposition, Eq. (|13|l . can be per- 
formed analytically. 

General symmetric states, = X^m=o a "i|^)j can be 
represented as 



*+j 7; 7 — I V Af-fe U/fe 

^N—kA^k.j 



N N-k k 

TN-kM^) = E a,m E E^ m > 

m=0 i=0 3=0 

(43) 

However, for these general symmetric states, the Schmidt 
coefficient matrix, c, is not sparse. 



E. Dynamics in the Symmetric Space 

One of the objectives of this paper is to treat dynam- 
ically generated entangled states, therefore this section 
briefly discusses the time evolution of symmetric states. 
It is straightforward to show that acting on a symmetric 
state with operators of the form, 



(44) 



preserves the exchange symmetry in the large basis, 
[o, Iljj] = 0, provided that the are identical. 

Using the symmetric state change of basis operator, 
Sn, elucidates the physical nature of the symmetric sub- 
space. For example, transforming any angular momen- 
tum operator of the form in Eq. 14411 to the small basis 
using Sn 



(45) 



produces the (N + 1) dimensional operator equivalent to 
the angular momentum for a single pseudo-spin J = N/2 
particle. This is because the symmetric subspace is com- 
posed of basis states \m) that correspond to the eigen- 
states of J z with J = N/2 (e.g., for two spins, the sym- 
metric subspace includes the triplet, but not the singlet). 

The dynamics of any symmetric state are confined to 
the symmetric subspace provided that the Hamiltonian 
can be expressed as a function of operators all of the 
form Eq. I|44l) . Given a symmetry-preserving Hamilto- 
nian, the dynamics can be completely simulated with the 
small symmetric basis. Explicitly, an infinitesimal step 
of evolution can be written 

|#(t + di)) = S N (l + iHdt)\V(t)) 

= S N \9(t)) + idtS N H\9(t)) 

= S N \*(t)) +idtS N HS N S N \*(t)) 

= \V(t))+idtH\%(t)) (46) 



where we have used |*(f)) = S N \^(t)}, H 
and |*(i)) 
symmetric) 



SnHS n , 

and = S N S N \^(t)) (because \&(t)) is assumed 
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For many experimentally motivated iV-particle spin- 
1/2 systems, it is possible to express states using the sym- 
metric subspace and the dynamics using only symmetry- 
preserving operators. The only time this efficient repre- 
sentation fails to apply is when the symmetry is broken 
or the system is divided (as we consider throughout the 
paper). For example, the spontaneous local decay of any 
one spin is sufficient to break the symmetry of Eq. (|16[) . 
Depending on the form of the decoherence, some sym- 
metry may be retained (e.g., the particle exchange sym- 
metry of Eq. UlSp). Other treatments have addressed the 
effect of such decoherence on parameters related to entan- 
glement, such as the degree of spin squeezing |ij IT5I Il6| . 

IV. ENTANGLEMENT PROPERTIES FOR 
REPRESENTATIVE SYMMETRIC STATES 

Given the large number of possible iV-spin states, even 
when restricted to the symmetric space, it is clear that a 
systematic, yet compact approach to characterizing mi- 
croscopic entanglement is necessary. Toward this end, 
we characterize a set of representative symmetric states 
with a limited combination of measures, including the re- 
duced state entropy, the entanglement of formation, and 
the logarithmic negativity. The families (described in de- 
tail below) that we have selected display diverse entan- 
glement behavior — they differ in their degree of entan- 
glement at different size scales and in their robustness to 
particle loss. Naturally, any set of representative states 
will be incomplete in some aspect; however, our goal is to 
provide a detailed picture of internal entanglement with- 
out an excessive number of representatives. 

In this section we address the relationship between the 
degree of entanglement and its robustness to particle loss. 
While it has been a longstanding conception that the 
most entangled states are simultaneously the most frag- 
ile, we demonstrate that this is not necessarily true. Un- 
der certain useful definitions of entanglement, it is pos- 
sible to find heavily entangled symmetric states that are 
simultaneously robust. Similarly, the most fragile states 
are not always the most entangled. We also demonstrate 
that restricting our analysis to the symmetric subspace 
does not preclude the potential for significant entangle- 
ment. 



A. Symmetric Reference States 

We now briefly describe several families of representa- 
tive symmetric states using the notation introduced in 
Section lull In addition to GHZ states and the W family, 
we introduce a new parametrized family termed "comb 
states" , which prove important in investigating the max- 
imal boundary of certain entanglement measures. 

Throughout the rest of this paper, all states are as- 
sumed to be symmetric. In the interest of simpler nota- 
tion, we will express symmetric states as p even when it 



is more efficient to compute entanglement measures us- 
ing their p representation. Tilde notation is used only for 
the |m) states. 



1. GHZ States 
The well known GHZ states [3^ can be written, 

|GHZ> = (|0> + \N))/s/2 (47) 

using the notation from Section IIIII The GHZ family 
is generally considered to be the standard example of a 
highly entangled state. In several different contexts it 
has become the common unit of entanglement currency. 
For example, as a particular Bell state, the GHZ state is 
the end goal of entanglement distillation protocols. 

However, the GHZ family fails to maximize a number 
of monotones, including the entanglements of distillation 
and formation for a given bipartite split. Unlike previ- 
ous treatments [l^, we choose to work with these mea- 
sures under which the GHZ is not a maximally entangled 
state. Certain other measures such as the N-tangle cor- 
rectly recognize the GHZ as containing the most true N- 
way entanglement 0, 0, EH , but our focus will remain 
on notions of strictly bipartite entanglement. Still, the 
most practical defining characteristic of the GHZ state 
is its fragility to particle loss; tracing out a single party 
destroys all of the internal entanglement. 



2. Dicke States 

An important family of states with completely different 
character is the set of symmetric states with integer m 
excitations (spins up), |m), where m = 0, ...,N. Of 
course, these states are also known as the Dicke states or 
the eigenstates of J z , where the notation | J, M) is used 
with J = N/2 and M = m - N/2 = -N/2, N/2. 
The W state [37j , which is defined as the symmetric state 
with one excitation, |W) = |1), is a particular member of 

this family. Notice that |m) and \N — m) have the same 
entanglement properties because one is equal to the other 
if the quantization axis is reversed. These states exhibit a 
high degree of entanglement for m = 1, . . . , N — 1, while 
the states of m = 0, TV are completely separable. The 
defining characteristic of the Dicke state entanglement 
is remarkable robustness to particle loss. It has been 
proven that |1) optimizes the concurrence when all but 
two spins have been removed |3S|. the extreme opposite 
of the fragile GHZ behavior. It has also been proven that 
for single copies the GHZ and W cannot be converted into 
each other with LOCC operations on the individual spins 
with certainty |37j . further emphasizing their difference. 
For additional discussions of the conversion properties of 
entangled states see [3^, 0, ^] . 
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3. Comb States 

A parametrized family of practical importance, which 
we call comb states, are defined as 

/27 N/s 

l C ( s )) = ViV £ W 2 + ™ s > (48) 

m=-N/s 

In the |m) basis, these states have a comb-like struc- 
ture with m-independent weighting for the non-zero el- 
ements which are spaced by s excitations. Because the 
comb states non-trivially explore the full support of the 
symmetric basis they may be expected to access regions 
of entanglement space where \fh) states are forbidden. 
Luckily, particular comb states with an optimized spac- 
ing s turn out to contain near maximal entanglement 
for bipartite splits of any symmetric ensemble as will be 
shown numerically and proven in the appendix. 
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4- Random States 

Another way to numerically explore the full symmet- 
ric space is randomly. We define a randomly generated 
state |R) = J2 m r m\fn} where the coefficients r m are com- 
plex Gaussian random variables with averages E[r m ] = 0, 
E[r m r„] = 0, and E[r^r„] = S mn /(N + 1). Note that 
this distribution of states is independent of the basis in 
terms of which we have chosen to define the random co- 
efficients r m . If we write |R) = J2 m r 'm(U\m)) in a new 
basis U\rh), where U is an arbitrary unitary, the new co- 
efficients r' m have exactly the same Gaussian distribution 
as the coefficients r m . As a result this distribution de- 
termines a measure on (unnormalized) vectors in Cjv+i 
that is invariant under unitary transformations. More- 
over E[(R|R)] = 1 so the states are on average normal- 
ized and in fact the distribution of norms becomes very 
sharply peaked around one as N — > oo. In this limit we 
can regard the states |R) as being drawn from the natural 
unbiased distribution of pure states. In practice we ran- 
domly select these vectors for a fixed finite N of interest 
and normalize. 



B. Pure State Entropy of Entanglement 

For an initially pure, iV-particle symmetric state, there 
are ^TV/ 2J possible ways to partition the system into 
two parts. With symmetric states we can replace the 
labelling of a particular partition {A, B} with the num- 
ber of spins in each partition {Na = k, Nb = N — k} 
where k = 1, . . . , [N/2\ . The entropy can then be com- 
puted from either of the reduced density matrices [23[ . 
PN-k = T^kPN or p k = Tr N ^ k p N : 

E(\9),{k,N-k}) = S(Tr k \*){*\) 

= S(Tr N _ k | *)(*|) (49) 



FIG. 1: Entropy of entanglement for representative symmetric 
states (described in Section II with N = 50 particles as a 
function of the dimension of the bipartite split, {k,N — k}, 
where k = 1, . . . , [N/2\. The unobtainable bound log 2 (fc + 1) 
is the entropy that could be achieved by a non-symmetric 
product of the two symmetric subsystems, {A,B}. Several 
representative states nearly achieve this maximum. 



It can be proved that the entropy is a monotonically in- 
creasing, concave down function of k in this range |42j . 
(From this point on, we implicitly assume the rounding 
of non-integer numbers such that |_A/2J is implied by 
N/2 and { [A/2J , [A72] } is implied by {N/2, N/2}.) 

In Section flll Al we emphasized that a symmetric state 
with N particles can be represented on the product space 
of two symmetric spaces with N — k and k particles 
(Sjy-fc ®Sfc). For all states within this space, the state of 
Eq. ( witn d A = k + 1) maximizes E(\$f), {k, N - k}) 
at log 2 (fc+ 1). However, this state is not symmetric with 
respect to the exchange of any two particles across the 
split. We are interested in finding the upper bound for 
the states in the space Sat, which are only a subset of 
states in E>N-k ® Sfe. It has been proven that the addi- 
tional restriction of overall symmetry constrains the max- 
imal entropy to be strictly less than log 2 (£;+ 1) except for 
N = 2,3,4 and 6 where states can be found that achieve 
this bound |l9j. Consequently, we refer to the bound 
log 2 (/c + 1) as the unobtainable bound for any k. 

Figure ^ shows a plot of E(\^>) , {k, N — k}) for several 
reference states and N = 50. Despite the fact that all 
states are forbidden from achieving the value log 2 (fc + I), 
some states come close to achieving this unobtainable 
bound. These include most randomly generated states 
and the comb states with s = \/2N. This naturally leads 
us to the question: what exactly is the minimum upper 
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bound for the split entropy of symmetric states and what 
states achieve that bound? 



1. Maximizing the Even Split Entropy 

Because the entropy is maximized by the most even 
split (k — N/2), we henceforth consider only this parti- 
tion. From the above discussion, we know that for N > 7, 
the entropy obeys the inequalities, 

E(\V), {N/2, N/2}) < E max {N) < log 2 (A/2 + 1) (50) 

Analytically locating the minimum upper bound 
Emax(N) (or the states that achieve it) is difficult, but 
a simple numeric approach turns out to shed some light 
on what we can expect. Figure shows the entropy 
of the even split entropy as a function of N for several 
families of states with the unobtainable upper bound for 
reference. Most families of states do not keep up with 
the scaling of this upper bound. 

For example, if N 3> m the states |m) (with m = 
1, . . . , N/2) can be shown to have entropies of 



E(\m),{N/2,N/2}) 



log 2 (m) 



(51) 



We also see that the largest of these scales as 
E(\N/2),{N/2,N/2}) » log 2 (iV)/2. 

Because of the factor of two, none of these states 
keep up with the Sjv/2 <8> Sjv/2 bound. However, if we 
explore the simplest possible states accessing more of 
the symmetric Hilbert space, we find something quite 
different. For large N (up to 600), the average en- 
tropies of random states, for example, numerically scale 
as « log 2 (7V/2 + 1) — 0.6. This indicates the remarkable 
fact that the symmetry constraint on the overall state 
does not limit the scaling of the maximal bipartite en- 
tanglement compared to that of the more general space 

Sat/2 <8> §AT/2- 

The comb states, optimized over the spacing s, are 
even more entangled. Numerically, we find that (for N 
up to 600) their entropies scale as w log 2 (-/V/2 + 1) — 0.3, 
when s ~ y/2N. Encouraged by this evidence, we were 
able to prove in the asympto tic l imit of large N that 
this family of comb states \C(s/2N)) does indeed scale as 
log 2 (N/2 + 1) — 8 where 8 is a constant of order unity (see 
appendix) . A similar proof for the random state scaling 
is probably possible. The fact that random states, and 
the optimized comb state, seem to nearly maximize the 
{N/2, N/2} entropy indicates that the set of states which 
scale similarly is of non-zero measure (i.e., this behavior 
is not atypical). 

Still we have not located the value of the true minimum 
upper bound and the form of the states that achieve this 
bound. Given the above results we expect it to have a 
similar scaling with a minimal offset, 8, for large N. 



C. Entanglement of Formation: Extremal Splits 

For any bipartite entanglement measure, we can con- 
struct even more possible splits if we choose (or are 
forced) to ignore some of the particles. Suppose we start 
with a symmetric state of N spins |$) and trace out 
spins until only N r remain. In this case, the new state 
PN r = TrAr_jv r (| v I')(^'|) will be mixed but still symmet- 
ric. We then have the possible bipartite splits {k, N r — k} 
with k = N r /2 [53|. 

For pure states, the entropy of entanglement for any 
bipartite split is equal to both the entanglement of forma- 
tion and distillation. Unfortunately, numerically calcu- 
lating either of these monotones is much more difficult if 
given an initially mixed density matrix. For negativities, 
we showed in section lill Cl that we can numerically calcu- 
late all bipartite splits {k, N r — k} for symmetric states, 
and we will demonstrate this ability in section llVDI For 
now we would like to deal with the extreme case of all 
but two spins removed (N r = 2). In section III Bl we 
stated that the entropy of formation Ef(\^), {1j 1}) is 
easily calculated for two spin mixed states through the 
concurrence. By discussing the relationship of the pair 
{E F (\^),{l,l}),E F (\y),{N/2,N/2})}, we can start to 
get a sense of the allowed relationship of entanglement 
across the extremes of size scales. We will refer to the 
splits {1, 1} and {N/2, N/2} as the extremal splits. 

Figure 03 displays Ep (\^), {1, 1}) for several refer- 
ence states. It has been proven that the W state |1) 
maximizes the concurrence, hence also the entanglement 
of formation (for all symmetric states) with a value of 
C(|l» = 2/N @||. Wang and M0lmer OH have 
shown that by using a similar formalism where the two 
spin concurrences are calculated from the moments of the 
entire state, analytic expressions can be derived for the 
concurrences of several families of symmetric states. In 
particular, for the Dicke states \m), and M = m — N/2, 
the concurrence is 

C{\M + N/2)) = 77 ^ 771 i — —{N 2 - 4M 2 



2N(N-1) 
- ^J{N 2 - 4M 2 )[(7V - 2) 2 - 4M 2 ]} 



(52) 



which gives the above result for the |1) state and also 
C(\N/2)) = l/(N - 1). In the large N limit these con- 
currences lead to the entanglements of formation 



M|i),{M}) 

E F (\Nj2),{l,l}) 



2 log 2 N + log 2 e 
N 2 

21og 2 2(A-l)+log 2 e 

(2(N-i)r 



(53) 
(54) 



The 1/N 2 scaling is due to the fact that the two spin 
state is constrained to be reduced from a larger sym- 
metric N spin state. In effect, one spin can only be so 
entangled with another when it is constrained to have 
the same relationship with all other spins. 

For many states Ef(\^), {1, 1}) is simply zero. The 
GHZ state, the comb state, and practically all random 
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FIG. 2: (A) Plot of the even split entropy of entanglement, E(p, {N/2, N/2}), for representative states as a function of the 
number of particles, N (which is also equal to the entanglement of formation and distillation). Note that the average entropy 
of 25 random states, \R), as well as the entropy of \C(V2N)} , nearly attain the unobtainable bound log 2 ([N/2] + 1). (B) A 
plot of the two-particle entanglement of formation, Ef{p, {1, 1}), as a function of the number of particles, N. The W state, 
jl), maximizes this entanglement measure, which quantifies robustness to particle loss. 



states have zero {1, 1} entanglement and do not contain 
the same degree of robust entanglement as the \fh) states. 
Furthermore the ordering of states shown in figure [2] A 
and B is reversed (with the exception of the GHZ). This 
leads to the question: what is the nature of the trade-off 
between the small and large scale entanglement of the 
extremal splits? 



Figure |3| shows each state as a point in the space of 
[E F (\*),{l,l}) t E F (\9),{N/2,N/2})] for N = 50. The 
line between |0) and |1) represents states which are a lin- 
ear combinationjrf these two states. The curve extend- 
ing from the \N/2) state to the vertical axis and up that 
axis to the comb state represents linear combinations of 
those two states. The forbidden regions of this space for 
symmetric states are unknown but we strongly suspect 
several properties of the boundaries. We conjecture that 
there are two regions where no states are allowed to ex- 
ist. First, in region (I), beneath the |0) «-> |1) line, no 
states are found, nor likely to exist. The reason for this 
is that to have any {1,1} entanglement there must ex- 
ist some degree of {N/2, N/2} entanglement. However, 
there must also be a region in the upper right (II) where 
no states exist. It appears that there exists a fundamen- 
tal trade-off between small and large scale entanglement: 
as the large scale entanglement of a state increases (and 
is above 1), the maximum allowable small scale entangle- 
ment will decrease. In addition, there is likely a critical 
value of the {N/2, N/2} entanglement above which the 
{1, 1} entanglement must be zero. 
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FIG. 3: Plot of accessible entangled states in a space 
that reflects the trade-off between the degree of entangle- 
ment (Ef(p, {N/2, N/2})) and its robustness to particle loss 
(Ef(p, {1, 1})). The degree of entanglement for the large split 
is also equal to the distillable entanglement. 



D. Negativities: Extension to All Splits 

Now that we have a better intuition for the relationship 
between the entanglement of the extremal splits, we can 
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more confidently approach the problem of understanding 
the large number of remaining splits. For N r spins re- 
maining, there are N r /2 splits of the form {k,N r — k} 
with k — 1, . . . , N r /2. If N r < N and the initial state is 
non-separable, the reduced state is mixed and one of the 
few computable entanglement measures available is the 
negativity. Even though it is a computable monotone, 
the negativity is not an entanglement measure with as 
much physical justification as the entanglement of for- 
mation or distillation. However, the logarithmic nega- 
tivity is an upper bound for the distillable entanglement 
|27| . With this in mind, we move forward and work with 
the logarithmic negativities as an indicator of potential 
entanglement. 
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1. Negativity of All Even Splits 

Before computing the negativities, we can use the 
properties of monotones to notice a few relationships be- 
tween the bipartite monotones of different splits. Tracing 
out a single spin is an operation that falls under LOCC, 
and any monotone X, including the negativity, can only 
decrease under such an operation, therefore 

E x {\*),{k-l,N r -k}) < E x (\^),{k,N r -k}) 
E x (\*),{k,N r -k-l}) < E x (\V),{k,N r -k}) 

For pure states, the most even split {N/2, N/2} gives 
the maximal entropy of entanglement [43] ■ We observe 
that this is also true for the most even splits of a reduced 
mixed state with N r particles remaining {N r /2, N r /2}. 
These observations motivate us to reduce the number of 
splits considered to only the even splits of a given N r . 
Figure H displays the quantity Em(\V) , {N r /2, N r /2}) as 
a function of N r for several reference states. The end- 
points of this plot give similar information about the ex- 
tremal splits as the previous description of entanglement 
of formation. Unlike the entanglement of formation, we 
can easily plot the intermediate splits for the logarithmic 
negativity. 

By the above inequalities, we know that each curve 
monotonically increases with N r . For reference we have 
included the plot of log 2 (AT r /2-|-l) which of course cannot 
be achieved, because each reduced state with N r spins 
remaining is constrained by the symmetry of the initial 
pure state. The space between this maximum and the 
space of all actual curves represents the entanglement 
'cost' of initial symmetrization. An unanswered ques- 
tion is for a given N r and N what pure state \$f) max- 
imizes £vv"(|*)){^V r /2,iV r /2})? What is this maximum 
as a function of N r and iV? These questions for both 
the negativities and other bipartite monotones are exten- 
sions of the problems encountered for the extremal splits. 
Again, we plot only the reference states and set aside the 
problem of fully characterizing the space of interest. 

First consider the GHZ state |0) + \N). As expected, 
this state is maximally fragile, starting at unity and drop- 
ping to zero as soon as one spin is removed. In direct 



FIG. 4: Plot of the even split negativity, 
Ej^{p, {N r /2, N r /2}), for representative symmetric states 
with N = 50, as a function of the number of particles 
remaining, N r , in the system. The inset plot highlights the 
particular robustness of the \N/2) Dicke state as measured 
by the negativity. This contrasts a similar analysis using the 
entanglement of formation, where |1) is most robust. 



contrast, the W state |1) starts at unity, but only slowly 
decays to zero as spins are removed and its logarithmic 
negativity remains finite for even N r — 2. The state 

\N/2) is, in some sense, an optimal trade-off between 
total entanglement and robust entanglement in that it 
starts reasonably high above unity at N r = N , but ap- 
pears to have maximal negativity below N r w N/2. The 
comb states (and random states), which have near max- 
imal total entanglement, are also a reasonable trade-off, 
especially compared to the extreme fragility of the GHZ 
state. 

The comb state and most random states still 'bottom 
out' with zero negativity (no negative eigenvalues of the 
partial transpose) below a critical N r . Because the log- 
arithmic negativity is an upper bound on the distillable 
entanglement, this must also be zero at these points. The 
size of this critical N r for a given state is another indi- 
cator of fragility of the entanglement (for the GHZ state 
it is the extreme N — 1). For the optimal comb state, 
the size of this critical value appears to scale only loga- 
rithmically in AT. Thus the comb states, despite having 
near maximal {N/2, N/2} entanglement, contain entan- 
glement that can withstand a huge amount of particle 
loss. 
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we can show 
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FIG. 5: (A) Plot of the inconsistent ordering of the reduced 
entanglement of formation, E F (p, {1, 1}), and the reduced 
logarithmic negativity, Ej^, {1, 1}), for Dicke states. The 
shaded region reflects the possible values for the distillable 
entanglement. (B) The large N scaling of the entanglement 
measures in (A) for the |1) (i.e., W) state. 



2. Ordering of Dicke States 

Given the fact that |1) optimizes the entanglement 
of formation of the {1,1} split, it may seem odd that 

\N/2) maximizes the negativity. Indeed, there is an or- 
dering issue here and the two monotones disagree on 
which of the reduced states is more entangled. See 
|22^ for a more complete discussion of ordering prob- 
lems with entropies and entanglement measures for two 
spin systems. Figure [S] displays the ordering problem be- 
tween EssQm), {1,1}) and E F (\m), {1, 1}) for N = 50. 
For N/2 >j^l, E F (Q), {1,1}) >J5 F Qj + l) t {l,l}) 
whereas E M (\j), {1, 1}) < Etf(\j + 1), {1, 1}), so the 
quantities are respectively decreasing and increas- 
ing with j. In fact the two curves will al- 
ways cross because, E F (\1),{1,1}) > E^Ql), {1, 1}) 
and E F (\N/2},{1,1})^ < E U {\N/2), {1, 1}). For 
large N 2 E F (\N/2), {1, 1}) a log 2 (iV)/2 < 

N 2 E A f(\N/2),{l,l}) ps iVlog 2 (e) where the approxima- 
tions can be shown both analytically and numerically. 

For N > m, N 2 Ej^(\fh), {1, 1}) flattens out to a con- 
stant as a function of N, while N 2 E F (\fh) , {1, 1}) con- 
tinues to grow logarithmically, as shown if figure |5j3 for 
|1). In this case, the entanglement of formation is signifi- 
cantly greater than the logarithmic negativity and hence 
also the distillable entanglement. So, for the state |1), 



E D (\1),{1,1}) < E„(\T), {1,1}) « w 

jp rh\ ji M\ ~ 2 log 2 N + log 2 e 
< E F {\1),{1,1}) ps — (55) 

All measures monotonically decrease with N, but the dis- 
tillable entanglement decreases at least logarithmically 
faster than the entanglement of formation. Similar state- 
ments are possible about any |m), with N 3> m. 

V. ENTANGLEMENT IN SYMMETRIC 
DYNAMICALLY GENERATED STATES 

Characterizing the reference states enabled us to quan- 
titatively identify the trade-off between the degree of en- 
tanglement and robustness of the state to particle loss. 
This relationship can be expressed as boundaries in the 
space expressed by the entanglements of formation for 
the extremal splits. With this relationship in hand, we 
are now able to address the question of where various 
dynamically generated states lie with respect to all ac- 
cessible symmetric states. 

For any given generation process, an important ques- 
tion involves exactly how entanglement forms within an 
ensemble 44]. In this section, we characterize spin- 
squeezed states, the most common experimental exam- 
ple of large scale entanglement. It has been shown that 
spin-squeezing (Eq. Q) is a sufficient condition for an 
iV-particle system to be entangled and the squeez- 
ing parameter also indicates in some sense the depth of 
entanglement [lj. ^ has also been demonstrated that 
spin-squeezed systems contain significant pairwise entan- 
glement [H,E3- However, little is known about the en- 
tanglement of squeezed states across all size scales or how 
they compare to the reference states from Section HVI De- 
scribing such states in terms of entanglement measures is 
intrinsically important, but also useful for understanding 
the more general class of symmetric entangled states. At 
the end of this section, we also briefly discuss the prob- 
lem of efficiently creating desirable states given speci- 
fied resources, allowable processes, and initially separable 
states. 



A. Spin Squeezed States 

The collective angular momentum operators of any 
multipartite spin state must satisfy the inequalities im- 
posed by their commutation relations. Let us assume 
without loss of generality that all subsequent states sat- 
isfy (J x ) = (J y ) = and {J 2 ) = mine((j|)) meaning x 
is the direction of the smallest variance perpendicular to 
the mean which points in the z direction. In this case, 
we use the uncertainty relationship 

(Jl)(j 2 y)> { -^- (56) 
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The characteristic feature of spin-squeezed states is that 
internal correlations between spins (i.e., entanglement) 
conspire to reduce the noise in one angular momentum 
component (x) at the expense of increasing the uncer- 
tainty in another (y). In particular, spin squeezed states 
satisfy the inequality, 



£2 = < , 



(Jz 



(57) 



States with a minimal squeezing parameter, £ 2 , are useful 
for reducing noise in many interferometric applications 
(e.g., atomic clocks). Using Eq. (j56|) and the fact that 
(J 2 ) < J 2 , one can show that, 



e> 



i 

N 



(58) 



where 1/iV is the Heisenberg limit. 
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1. Squeezing and Entanglement 

We choose to generate near optimally spin-squeezed 
states l^g) by applying the counter-twisting Hamilto- 
nian H ct = (J 2 . — Jz)/i to an initially polarized sample 
l^o) = |0) (with £ 2 = 1) for the length of time t/v needed 
to minimize £ 2 0- This process does not produce opti- 
mally squeezed states (see H2), but in the large N limit 
it creates states which very nearly obtain the minimal 
value of £ 2 . The time it takes to reach the minimum of 
<f oc 1/N for large N is t N 0.2 log 2 (A)/JV 0. Hence- 
forth, time is scaled such that the optimal spin squeezing 
time tj\r = 1. We will ignore the small difference be- 
tween the achieved and optimum spin squeezing so that 
we may examine the production of entanglement as the 
state evolves in the most simple way. Interestingly, an ef- 
fective counter-twisting Hamiltonian can be experimen- 
tally realized through the QND-detection and feedback 
rotation scheme of |8(. 

Figure shows this evolution for a state with N = 50 
spins. The x and y means remain zero for the entire 
evolution, while (J z ) decays from completely polarized 
through zero. For small numbers of spins, the state will 
quickly re-cohere and become completely polarized (sep- 
arable). For large numbers of spins, the dynamics be- 
come highly disordered after the mean decays through 
zero, and the re-coherence time grows much longer. Af- 
ter becoming maximally squeezed, the internal entangle- 
ment continues to grow, but the spin-squeezing rapidly 
gets worse because of the reduction in the mean. The 
entanglement of formation for the largest and smallest 
even splits [E F (p{t),{N/2,N/2}) and E F (p(t), {1,1})} 
are also shown normalized by their own initial local max- 
imum. 

The small scale entanglement {1,1} reaches its peak 
before the large scale entanglement {N/2, N/2} does. If 
we analyze the relative rate of growth of the different 
scales of entanglement at early times we see an intuitive 



FIG. 6: Spin squeezing evolution for a system of N = 50 spin- 
1/2 particles evolving by the counter- twisting Hamiltonian as 
measured by the squeezing parameter, £ 2 . The time is scaled 
such that maximal spin squeezing occurs at t — 1. The mean 
J z and the entanglements of formation are all independently 
normalized by their own maximum in the time period shown. 
Notice that the small scale correlations, Ef(p, {1, 1}), peak 
before the large scale correlations, Ef(p, {N/2, N/2}), as the 
squeezing evolves. 

ordering. Figure |7| shows the small-time logarithmic neg- 
ativities (for all even splits) and the entropy of forma- 
tion (for the extremal splits) normalized by their respec- 
tive maxima over that interval. As the state becomes 
squeezed, the {1,1} correlations form first, followed by 
the {1,2}, then the {2,2}, and so on up to {N/2, N/2}. 
This observation suggests that small scale correlations 
typically peak earlier than larger scale correlations when 
evolving under quadratic Hamiltonians. 

Another observation is that for small times, the state 
gets progressively more entangled in the sense of ma- 
jorization |45j |. In other words, the eigenvalues of 
Trk(p(t + dt)) are more disordered than the eigenvalues 
of Trfc (/?(*)) for all k < N/2 and small t. Thus, despite 
certain ordering difficulties with various entropies, the 
entanglement of any split is strictly increasing initially. 

It is also important to quantitatively compare the en- 
tanglement measures for spin-squeezed states and the 
symmetric reference states. Figure |SJ^ shows the even 
split entropy Ep(\^ 6 ),{N/2,N/2}) of the optimally 
squeezed state as a function of N. From a numerical 
fit, we find that 

E F (\Ve),{N/2,N/2}) w 0.461og 2 (A0 - log 2 (e) (59) 

For smaller scale entanglement, Figure [SJ3 displays the 
two-spin entropy E F (\^^) , {1,1}). The values approach 
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FIG. 7: Entanglement measures for a system of N = 50 
spin- 1/2 particles evolving under the influence of a counter- 
twisting spin-squeezing Hamiltonian. The time is scaled 
such that the squeezing parameter achieves its minimum at 
t = 1 (the small-time evolution is depicted) and all en- 
tanglement measures are independently normalized by their 
own maximum in the time period shown. The entropy 
of formation, Ef(p, {A, B}), is shown for the extremal bi- 
partite splits, {1,1} and {N/2, N/2}, while the logarith- 
mic negativity, Ej^f(p,{A,B}), is depicted for the partitions, 
{1, 1}, {1, 2},{2, 2},..., {N/2 - 1, N/2}, {N/2, N/2}. It can 
be seen that small-scale correlations tend to peak before 
their large-scale counterparts; the entanglement measures are 
strictly ordered according to the number of particles remain- 
ing N r . 



but never exceed the curve for \N/2). Indeed, it can be 
shown that in the large N limit, the two-spin concur- 
rence scales identically for the two states: C(|iV/2)) ss 
C(|^)) « 1/N, thus the entanglements of formation 
must also converge. 

For a specified number of particles (N = 50), the fam- 
ily of states generated by applying the counter-twisting 
Hamiltonian to a polarized sample are displayed in the 
[E F (p(t), {N/2, N/2}), E F (p(t), {1,1})} space of figure 
Again the small scale entanglement grows faster than 
the large scale entanglement, but eventually decays to 
zero as the large scale entanglement takes over. The 
disordered nature of the counter-twisting Hamiltonian 
dominates at long times as the value of the large scale 
entanglement diffuses and the small scale entanglement 
remains near zero. In contrast, the application of a twist- 
ing Hamiltonian H t — (which, unlike the counter- 
twisting Hamiltonian, creates squeezed states with a ro- 
tating axis of squeezing) is seen to be much more periodic. 
The states it generates are similar to the counter-twisting 
states initially, but they eventually converge to the GHZ 



state then return along the same trajectory. 

The entropies of extremal splits ({1, 1} and 
{N/2, N/2}) capture much of the character of a 
many-particle entangled state, but there are of course 
a large number of other bipartite splits to consider. 
The introduction of the information contained in all 
other splits potentially brings up more interesting 
entanglement characteristics. As in Figure 0] we can 
efficiently calculate all even split bipartite logarithmic 
negativities for large number states as they become 
spin squeezed. The characteristic of early small scale 
entanglement being transformed into subsequent large 
scale entanglement during the course of the evolution is 
again apparent. Nonetheless, for this particular case, 
the intermediate splits do not provide a considerable 
amount of additional insight compared to that from the 
extremal splits. 

2. Squeezing Under Particle Loss 

We now address how the spin-squeezing parameter be- 
haves under particle loss. Given the expectation values 
of a set of operators on a symmetric density matrix, it is 
simple to determine the moments of the same state with 
a certain number of particles removed. 11 pn is symmet- 
ric, so are all of its reduced density matrices pN r where 
1 < N r < N. Given single particle operators we know 
that 

Tr N (oi ■ ■ ■ Ojp N ) = Trjv_i(oi • • • o^pjv-i) (60) 

assuming the indices of the operators are not the ones 
traced out. With this observation and the fact that for 
symmetric states 

(J Z ) N = NTr N (j z ,ip N ) (61) 
(J z )n-i = (N - l)Trjv-i(j z ,iPjv-i) (62) 

we find 

N - 1 

(J.)n-i = ~^{J Z )n (63) 
Similarly, it is easy to show 

{JI)n-x = ^^{JI)n+\ (64) 

Using these relations and taking the large N limit we 
find that the spin squeezing parameter of a state with N r 
spins remaining ) is dependent on the initial squeez- 
ing parameter (£jy) and polarization of the state with all 
spins remaining in the following way 

& r =tf + (&-tf)^£y (65) 

where £ = N 2 /(4( J z )%). 

The inset of Figure [§] shows the spin-squeezing behav- 
ior for N = 300 as a function of time. Considering only 



15 




50 100 150 200 250 300 1 2 3 4 5x10^ 

Number of Particles (N) Ef(p,{1 , 1 }) 



FIG. 8: (A) A plot of the even split entanglement of formation (and entropy), Ef(p, {N/2, N/2}), for a system of N spin-1/2 
particles evolved under a counter-twisting spin-squeezing Hamiltonian. The state l^) minimizes the squeezing parameter, 
£ 2 . (B) A similar plot using the scaled entanglement of formation, N 2 Ef(p, {1, 1}), for a system with all but two particles 
removed. (C) The time evolution of states evolving under both the counter-twisting Hamiltonian (H c t = (J+ — J-)/i) and 
twisting Hamiltonian (Ht = ) in the space of extremal split entanglement. 



the time when the state of all the spins is maximally 
squeezed (t = 1), we plot the spin squeezing parameter 
as a function of number remaining in FigureEl which be- 
haves according to Eq. Ij65(l . For this finite number case, 
the spin-squeezing is lost after some fraction of the spins 
are removed. As N goes to infinity though, ( J z ) m — > N /2 
(£i - > 1) an d £jv —> so all spins need to be removed for 
the state to completely lose its spin squeezed character. 

In a similar analysis, Simon and Kempe fl5| have 
shown that spin-squeezed states remain squeezed until 
more than 29% of the particles have depolarized. Thus 
spin-squeezed states are robust to both par ticle loss and 
dephasing with constant N. (See 0, Il6j ] for a more 
complete treatment of how the spin-squeezing parameter 
behaves under continuous generation and decoherence.) 
However, robustness to particle loss and dephasing do not 
necessarily imply each other because the GHZ states are 
remarkably robust to local depolarization [l]| , but obvi- 
ously maximally fragile under particle loss. The complete 
relationship between robustness to particle loss and de- 
phasing is an interesting direction for further research. 



B. Generating Entangled States 

Instead of characterizing what states a particular pro- 
cess produces, consider the reverse problem of determin- 
ing the process necessary to generate a desired state from 
an initially separable state. The completely polarized ini- 
tial state (|0) or \N)) is usually chosen both because it is 
completely separable at all levels and it is easily prepared 
in the lab (e.g. via optical pumping). 



It can be shown that given such a state and access to 
Hamiltonians of the form J x , J y , J z , and J^, i.e., the gen- 
erators are all rotations plus a single non-linearity, one 
can produce any symmetric state by an algorithm that 
switches between the Hamiltonians in time [46) . Unfortu- 
nately, proving this statement does not necessarily lead to 
the most efficient way to create a particular state. Know- 
ing which states are prohibitively expensive to produce 
is an important experimental question. An interesting, 
but difficult, way to characterize a state is by quanti- 
fying the resources needed to create that state given a 
certain set of generators. For example, one could define 
a cost metric which is a function of how many times the 
Hamiltonians must be switched and the length of time 
necessary to produce a particular state. 

Of course, all of these issues are context specific, but 
we can summarize certain results. Simply observing 
what the application of a particular Hamiltonian pro- 
duces is a first step. The counter-twisting Hamiltonian 
presented earlier produces optimal squeezing but docs 
not produce any recognizable reference state (since the 
dynamics for large N becomes highly disordered for long 
times). A one- axis twisting Hamiltonian of produces 
some squeezing which does not scale optimally 2] . How- 
ever, the time dependence of the entanglement produced 
by this Hamiltonian is much more periodic and ordered 
than the counter- twisting version. In fact, it produces 
the GHZ state halfway through its period as is indicated 
in figure rvfolmer and S0rensen have proposed a 
robust scheme for generating the GHZ state of many hot 
ions taking advantage of this effect. 

Unanyan et al. have shown that by using adiabatic pas- 
sage and energy level navigation methods one can pro- 
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FIG. 9: Spin squeezing reduction due to particle loss for a 
system initially containing N = 300 spin- 1/2 particles. The 
lower right corner of the plot reflects the minimum value of 
the squeezing parameter, £ 2 , achieved via a counter- twisting 
Hamiltonian. As particles are removed from the optimally 
squeezed system, the value of £% r moves up and left along 
the plotted line, eventually crossing unity for finite number 
systems. The inset shows the time evolution of £ 2 for different 
numbers of remaining particles, N r . 



duce the GHZ state and all |m) states |48|. However, 
it remains unclear what the most efficient method is to 
generate these states, or the bipartite entropy maximiz- 
ing states presented here, in the asymptotic limit of large 
N. 



VI. CONCLUSION 

In this paper we analyzed the microscopic structure 
of entanglement and its robustness to particle-loss for 
many-particle symmetric states. Our approach pro- 
ceeded by comparing the features of dynamically gen- 
erated squeezed states to a collection of symmetric rep- 
resentative states, including the GHZ and Dicke states, 
as well as random states and a new family that we de- 
fine. In order to perform the analysis, we selected sev- 
eral bipartite entanglement measures: the reduced en- 
tropy of entanglement, the entanglement of formation, 
and the logarithmic negativity. By computing these bi- 
partite measures for all possible reductions and partitions 
of the systems, we were able to construct a picture of 
multi-scale entanglement. 

Our analysis benefitted from simulations of many- 
particle systems. The computational results helped to 
bolster physical insight and provide a starting point for 
analytically treating the asymptotic scaling of various 



entanglement measures. In order to circumvent the ex- 
ponential scaling of the density matrix for arbitrary N 
particle states, we restricted our analysis to the symmet- 
ric subspace. In Section IlIII we developed machinery for 
computing the above entanglement monotones for sym- 
metric states in a computationally efficient manner. As a 
result, our simulations were capable of handling systems 
with N ~ 10 3 particles without making any dynamical 
approximations. 

In Section IIVI we characterized the entanglement of 
the representative states in detail, focusing on the trade- 
off between those states that maximize the entanglement 
measures and those that are robust under particle loss. 
We also analyzed several important ordering issues be- 
tween the different measures. A key point we stress is 
that fragility is not necessarily a property of highly en- 
tangled states. With the analysis in Section IIVI it was 
possible to address the evolution of microscopic entan- 
glement in dynamically generated spin-squeezed states. 
Hopefully this work helps clarify the otherwise vague 
statement that "spin-squeezed states are massively en- 
tangled." 

From this work, we anticipate several future direc- 
tions. First we plan to consider less restrictive sym- 
metry classes, particularly the symmetry of Eq. 1)15(1. 
This symmetry is preserved during the unconditional 
evolution of an ensemble under a uniform symmetry- 
preserving Hamiltonian and local dephasing. For cer- 
tain cases where the emission from the particles does not 
physically distinguish different particles, the symmetry 
may also be preserved under conditional evolution. In 
order to perform such an analysis, it will be necessary to 
exploit the commutant algebra and representation theory 
of the permutation group jl^ ■ Preliminary investigation 
suggests that it will be possible to treat the full permuta- 
tion group in a manner that scales polynomially, rather 
than exponentially, with the number of particles. 

A more straightforward goal is to generalize the treat- 
ment of this paper to particles with more than two levels. 
For example, we would like to describe the entanglement 
within an ensemble of Cesium atoms, where each atom 
can occupy the nine magnetic sub-levels of the F = 4 
ground state. 

Regarding dynamically generated states, it is possi- 
ble to further simplify the description of entanglement at 
small times. For any initially polarized state experiencing 
a quadratic Hamiltonian, the state and relevant entangle- 
ment measures can be described in terms of the Gaussian 
moments alone for short times. Deriving the exact form 
of this low-dimensional parametrization of entanglement 
measures is of experimental interest. 

Finally, an important challenge is to develop tech- 
niques for efficiently generating the reference states dis- 
cussed in this paper, including those with near maximal 
entanglement, such as the comb states. Here, we hope to 
stress that theoretical treatments of many-particle spin 
systems are most beneficial when they adopt methods 
that can be experimentally implemented. 
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APPENDIX: SYMMETRIC STATE ENTROPY 
SCALING 

Proposition 1. There exists a lower bound for the max- 
imum attainable symmetric state entropy that asymp- 
totically scales as the maximum entropy for states in 

e S N ,5>0,N* >0:VN>N*, (A.l) 
log 2 (N/2+l)-E(\V),{A,B})<5 

Proof. The proof proceeds by constructing a sym- 
metric state whose even split reduced entropy, 
E(\^>) , {N/2, N/2}) can be expressed as the asymptotic 
series, \og 2 (N/2 + 1) + S. In order to obtain this series, 
we express the entropy in terms of the Schmidt coeffi- 
cients from Eq. I|42|l by employing Stirling's formula. 
Computing the residuals that are incurred by Stirling's 
approximation, we obtain a bound for 6 and demonstrate 
that it is asymptotically constant in N, i.e., 5 ~ 0(1). 

Consider the family of |C(s)) states, defined in Section 
IIV A 31 whose Schmidt decomposition, according to Eq. 
l|4"2jl. is given by, 



N/s 



\C(s)) = A 



N-k k n X 
\ - \ - ^N,mON_ +ms i+j 



m=-N/s i=0 j=0 



CjV-fe,iCfc, 



i)N-k\j)l 



(A.2) 

where A = y/2s/N. We wish to choose the value of s in 
|C(s)) such that the matrix, c, becomes block-diagonal 
in the large ./V limit, which will happen provided that the 
\fh) contributing to \C{s)) are sufficiently separated. For 
an orthogonal Schmidt matrix, Eq. I|13|) can be solved in 
closed form and the total entropy is a weighted sum of 
the entropies contributed by each participating, \fh). 

The required separation between nonzero \fh) in \C(s)) 
as well as their contributing entropies, can be found by 
considering c m (i), the c%j matrix elements as a function 
of i for a given value of m. This leads to the distribution, 



Cm (^) 



1 



N-k 



(A.3) 



which can be approximated using Stirling's formula, 

1 



log n\ = n log n — n + 
to obtain the expression. 

c 2 m (i) = 2-<™+*)ex P 



27m + — + 0(n- 2 ) (A.4) 
lZn 



i 2 — i + m + to 2 



12i 2 m — 12im 2 



x (A.5) 
+ 0(m- 2 ) 



for c m (i) as N — > oo. The residual terms are of order 
m . 

Unfortunately, it is not known how to evaluate the 
entropy of Eq. (|A.5J) because the discrete sum, 



S(m) = - ^ c 2 m {i) log 2 c 2 m {i) 



(A.6) 



i=0 



cannot be expressed in closed form. However, the mo- 
ments of c m (i) can be computed analytically, 



i = (i = 



cr 2 = (i 2 ) - (i) 2 = —m(m + 2) 



(A.7) 
(A.8) 



and it can be seen that all higher cumulants vanish in 
the large N limit. 

The entropy contributed by each |m) in |C(s)) when 
m is large can be computed by approximating the c^(z) 
as Gaussian, 



exp 



-1 V3(i-m/2) 2 
4m yjm(m + 2) 



0(m 



and transforming the summation in Eq. (|A.6Jl into an 
integral, 



S m -e=- J c m (i) log 2 cf m {i) di +1 0(m )di 

(A.10) 

where the error term, e, must be obtained using the 
Abel-Plana procedure [5f| for computing the difference 
between a discrete sum and its corresponding inte- 
gral. An upper bound for the integral over the resid- 
ual 0(m~ 2 ) terms can be found to converge to, e = 
\/27rm~ 2 exp(— m 2 )/128, by considering the asymptotic 
series of the T function |5lT |. 

The resulting entropy, in the large m limit, with the 
Abel-Plana corrections and the error bounds from higher 
order Stirling terms can be shown to have the form, 



(All) 



e -l/4m m l/4( 2 + m )3/4 

m ~ 3 1 /4 v /2m(2 + m) log 2 

+^ 3 (0, m-i e- 1 / 4 ™) - it?4(0, 1/0F) + e 

where $„(«, v) is the elliptic theta function of order n 
[5l| . Computing the entropy of the full state requires 
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evaluating the discr ete s um over to. We chose the comb 
spacing to be s = \2N based on Eq. I|A.8(I . Performing 
the final sum leads to, 



1 



E(\V), {N/2, N/2}) = (1+— )=) log 2 N+0(Ne- VW )+e 



2VN 

(A.12) 

where the first residual term reflects the finite overlap 
of the c^(i), i.e. corrections that arise because c is not 



perfectly block-diagonal. 

Finally it is possible to express the reduced entropy as 
an asymptotic series, 

E^), {N/2, N/2}) -> log 2 (iV/2 + 1) - 0(1) (A.13) 

which completes the proof. 

□ 
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We work with the convention that for any bipartite split 
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of a system, {A,B}, the symbol A is equivalent to a 
list of the particles residing on its side of the split, and 
likewise for B. Na and Nb are the number of particles 
on each side of the split; (La and ds are the Hilbert space 
dimensions of each side, generally 2 Na and 2 Nb for non- 
symmetric collections of spin- 1/2 particles. 



[53] We use the convention that for any bipartite measure X if 
j + k is less than the number of spins N in |^}, the possi- 
bly mixed state that the measure Ex(\^), {j, k}) should 
act on is actually Trj V -j-fe|'l/}('I'|. For the concurrence 
we use the similar convention C(|*}) = C (Tr W -2 ] *> {* | ) • 



